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SOME RECENT CONTRIBUTIONS TO THE STUDY OF 
TRANSITION AND TURBULENT BOUNDARY LAYERS 1 
By Hugh L. Dryden 


SUMMARY 


The first part of this paper reviews the present Btate of the prob- 
lem of the instability of laminar boundary layers which has formed an 
important part of the general lectures .by von Karman at the first and 
fourth Congresses and by Taylor at the fifth Congress. This problem may 
now be considered as essentially solved as the result of work completed 
since 1938. When the velocity fluctuations of the free-stream flow are 
lees than 0.1 percent of the mean speed, instability occurs as described 
by the well-known Tollmien-Schlichtiiig theory. The Tollmien-Schlichting 
waves were first observed experimentally by Schubauer and Skramstad in 
1940. They devised methods of introducing controlled small disturbances 
and obtained measured values of frequency, damping, and wave length at 
various Reynolds numbers which agreed well with the theoretical results. 
Their experimental results, were confirmed by Liepmann, Much theoretical 
work was done in Germany in extending the Tollmien-Schlichting ; theory to 
other boundary conditions, in particular to flow along, a porous wall to 
which suction is applied for removing part of the boundary layer.' 

The second part of this paper summarizes the present state of knowl— . 
edge of the mechanics of turbulent boundary layers, and of the methods 
now being used for fundamental studies of the turbulent fluctuations in 
turbulent boundary layers. A brief review is given of the semi-empirical 
method of approach as developed by Buri, Gruschwitz, Fediaevsky, and 
Kalikhman. In recent years the National Advisory Committee for Aeronautics 
has sponsored a detailed study at the National Bureau of Standards of tte 
turbulent fluctuations in a turbulent boundary layer under an adverse 
pressure gradient sufficient to produce separation. The aims of this in- 
vestigation and its present status are described. 


1 Paper presented at the Sixth International Congress for Applied 
Mechanics, Paris, September 22-29, 1946. 
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INTRODUCTION 


Since the last meeting of the International Congress for Applied 
Mechanics in 1938, there have been important advances in our understand- 
ing of some of the most fundamental aspects of fluid mechanics, in 
spite of the fact that the major attention of scientists was devoted to 
war research during most of the period. The results of this basic re- 
search have been available to a limited number of workers in confidential 
or secret reports. These reports have now been declassified for the 
most part and the exchange of information between scientific workers is 
increasing as evidenced by the convening of this Congress. Nevertheless, 
it would be presumptuous of me to assume that I am familiar with all the 
scientific work in fluid mechanics in all countries, even in the rela- 
tively narrow field which forms the subject of this lecture. Apologies 
are therefore made to any colleague who may find his accomplishments 
overlooked. Of necessity a large part of the paper describes work done 
in my own country and in particular at my om laboratory with the coop- 
eration and financial assistance of the National Advisory Committee for 
Aeronautics. Acknowledgment is made to G. B. Schubauer, H. K. Slcramstad, 
P. S. ICLebanoff, and W. Squire, who carried out the experimental work at 
the National Bureau of Standards, and to Br. George W. Lewis, Director 
of Research of the National Advisory Committee for Aeronautics, who gave 
permission to use some data not hitherto published. 

This paper reviews the present state of knowledge of the mechanics 
of boundary— layer flow. As is well known, L. Prandtl (reference l) in 
1904 introduced the concept of boundary layer, a region of small thick- 
nesB near the surface of an object immersed in a fluid stream, or moving 
through a fluid, within which the speed of the fluid relative to the 
surface rises within a comparatively short distance from zero at the 
surface to a value comparable with the relative speed of the body and 
the fluid at a great distance. The rate of change of speed with dis- 
tance within the layer in the direction of the normal to the surface is 
at least an order of magnitude higher than in the flow outside the layer. 
Beginning about 25 years ago it has been possible to study the distribu- 
tion of speed and the character of motion within the boundary layer. 

It was soon realized that the state of flow in the boundary layer 
may be either laminar or turbulent corresponding to the two states of 
flow in a pipe which had been known since the time of Osborne Reynolds. 
Every layman has observed these two states of flow in the smoke rising 



NACA TN Mo. 1168 


3 


from a cigarette in a quiet room, that is, the smooth flow in layers, with 
lazy oscillation of the stream as a whole near the cigarette and the sud- 
den transition to a confused series of whirls and eddies at same distance 
above it. The essential distinction "between these two types of flow has 
nevertheless been the subject of long and profound discussion at many 
of our previous meetings. 

The key to many problems in the mechanics of fluids is to be found 
in the state of motion in the boundary layer. In a lecture at Toulouse 
in June 1943 M, J. Kampe de Far let (reference 2) reviewed the signifi- 
cance of boundary layer studies to practical aeronautics under the title 
"Un probleme-cle de l'aeronautique: 1 ! etude de la couche limite," giving 
special attention to the reduction in the drag of an airplane wing by 
smoothing the surface and maintaining laminar flow. The whole question 
of scale effect, the effect of wind tunnel turbulence on wind tunnel 
measurements of aerodynamic forces and moments on aircraft models and 
on control surfaces of aircraft models, flow separation, and wani mum 
lift coefficients — all these and many Bimilar questions can be ade- 
quately understood only in terms of knowledge of the state of flow in 
the boundary layer. 

Even at speeds approaching the speed of sound, the recent work of 
Ackeret (reference 3) and Liepmann (reference 4) has shown that the flow 
pattern, including the location and form of shock waves, depends markedly 
on the state of flow in the boundary layer. There is thus a Scale effect 
even at high Mach numbers which can be understood only in the light of the 
state of flow in the boundary layer. 

It is however not the intention in this paper to pursue in detail 
the effects of boundary layer flow on the general flow pattern and the 
consequent effects on aerodynamic forces and moments e We. shall restrict 
our attention to the boundary layer itself with a view to understanding 
the phenomena occurring within it. 

The equations of laminar flow in a boundary layer were given in the 
original paper of Prandtl (reference 1 ) in 1904 and a solution for the 
special case of flow along a thin flat plate parallel to the direction 
of flow and without pressure gradient was given by Blasius (reference 5) 
in 1908. Experimental data on the velocity distribution, were obtained 
by Burgers and Van der Hegge Zijnen (reference 6 } = in 1929 and later by 
others which confirmed the theoretical computations. It is certain 
that the equations are correct and hence that the fundamental mechani— ■ 
cal principles are fully understood. For that reason no attempt is 


\ 



k 


MCA TN No. 1168 


made here to review the many papers relating to laminar hound ary layers. 
Reference is made to the hook "Modern Developments in Fluid Dynamics" 
hy S. Goldstein and others (Oxford, 1938) especially to chapter IV, and 
to the papers listed in the bibliography under the heading, T-nmir^-r 
Boundary Layer . 

Thus more specifically this paper deals with two problems: (l) 
Stability of Laminar Flow in a Boundary. Layer, and (2) Mechanics of the 
Turbulent Boundary Layer. The general plan of procedure will be first 
to review briefly the status of each problem at the time of the last 
Congress in 1938. and then to present the important developments since 
that time . 


STABILITY OF LAMINAR FLOW IN. A BOUNDARY LAYER 
Status at 1938 Congress 


Cause of transiti on.- The problem of transition to turbulence in 
boundary layers received much attention at the last Congress. G. I, 
Taylor devoted nearly half of his general lecture to the subject and 
four additional papers had some bearing on the problem. There were at 
that time two schools of thought. On© school thought that transition is 
the result of a definite instability of the laminar boundary layer in 
which infinitesimal disturbances grow exponentially. The other thought 
that the laminar flow is stable for infinitesimal disturbances but that 
transition occurs as a result of external disturbances of finite 
magnitude. 

Many investigators have attempted to solve the theoretical problem 
of the stability of laminar flow by determining what conditions are nec- 
essary to cause small disturbances in the form of velocity variations to 
increase with time. This work dates from that of Lord Rayleigh (refer- 
ence 7) in i860. About 17 years ago the specific problem of the stabili- 
ty of the laminar flow in a boundary layer near a thin flat plate without 
pressure gradient was studied by W. Tollmien (reference 8) at Gottingen. 
The problem was idealized by assuming a layer of constant thickness ex- 
tending to an infinite distance in both directions within which the dis- 
tribution of mean velocity was that computed by Blasius , According to 
the results of Tollmien^ computation small disturbances in velocity 
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of vave length lying in a certain region would "be amplified whereas dis- 
turbances of shorter or longer wave length would be damped, provided the 
Reynolds number of the boundary layer was greater than a certain value. 

The calculations were repeated and extended by H. Schlichting (references 9 
and 10) in 1933 and 1935. The amplified disturbance was assumed to grow 
until it caused a breakdown of the laminar flow. These were the views 
of the theoretical workers of the Gottingen school in 1938. 

The experimental workers could find no evidence in 1938 to support 
the Tollmien-Schlichting theory. In 1935 my colleagues and I had stud- 
ied the flow near a flat plate experimentally in continuation of the 
work described at the 193^ Congress for Applied Mechanics and failed to 
find any evidence of the development of the amplified disturbances de- 
scribed by Tollmien and Schlichting. Experiments by the groups at the 
California Institute of Technology, Massachusetts Institute of Technol- 
ogy, and Cambridge University, England, likewise failed to Bhow the 
amplified disturbances. Hence the views of these groups were that the 
experimental evidence was definitely against the instability theory. 


Effect of free-stream turbulence .— An alternative view was that 
the. presence of finite disturbances in the free stream was the principal 
factor producing transition. . A theory was given at the 1938 Congress by 
G. I. Taylor of the effect of free-stream turbulence . He assumed that 
transition occurred as the result of momentary separation in the regions 
of adverse pressure gradient associated with the turbulent fluctuations 


at sufficiently large values of the Karman— Pohlhausen parameter S 2 /(iU x 
dp/d x, where 5 is the thickness of the boundary layer, u the vis- 
cosity of the fluid, U the speed of the fluid in the free stream, and 
dp/dx the pressure gradient. On this assumption Taylor found that the 


Reynolds number at transition was a function of 


uVDN 
U \L/ 


D y /s 


where is 


the intensity of the turbulence, L is its scale, and D is a reference 
dimension of the body. This relation had been experimentally verified in 
the Cambridge laboratory for flat plates, and in our laboratory for spheres 
and for an elliptic cylinder. 

L. Schiller presented at the 1938 Congress a study of transition in 
a pipe, attributing transition to the eddies originating at the entrance 
of the. pipe and giving a quantitative theory. 


Effect of curvature .- In addition to the experiments of G. I. Taylor 
(references 11 and 12) on the flow between rotating cylinders confirming 
the greater stability of flow over a convex as compared with a concave 
surface, the results of Clauser (reference 13) on the flow past a curved 
plate were available. These agreed qualitatively with the rotating 
cylinder result. but there were some quantitative differences. 
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Effect of pressure gradient .— The discussion of the effect of pres- 
sure gradient at the 1938 Congress dealt mainly with the possible de- 
stabilizing effect of .adverse pressure gradients. The information then 
available indicated that the boundary layer might be non-turbulent until 
separation occurs and that turbulence set in at this point at high 
Reynolds numbers and at some distance downstream from it at low Reynolds 
numbers . 

.fluc tuations in the laminar boundary layer .- It was suggested that 

the slow fluctuations in the laminar boundary layer reported at the 1934 
Congress were probably fluctuations in thickness of the layer. Since the 
components of the fluctuations are not correlated, there are no shear 
stresses, and the distribution of mean velocity is unaffected by the 
presence of the fluctuations. 


Discovery of Tollmien-Schlichting Oscillations 

In 1940 a research program was undertaken at the National Bureau of 
Standards with the cooperation and financial assistance of the National 
Advisory Committee for Aeronautics to investigate the effectiveness of 
damping screens in reducing wind-tunnel turbulence and to study transi- 
tion on a flat plate down to the lowest attainable turbulence level. 
Turbulence levels as low as 0.02 percent = 0.0002) were obtained. 

The Reynolds numbers for transition increased steadily as the turbulence 
was reduced to 0.08 percent and thereafter remained constant, as shown 
in figure 1. 

In the course of this work it was decided to repeat the earlier 
observations of fluctuations in the laminar layer under conditions of 
low turbulence and in August 1940 my colleagues, G. B. Schubauer and 
H. K. Skramstad, obtained the records of the velocity fluctuations in 
the boundary layer of a flat plate by means of a hot wire~anemometer 
as shown in figure 2. The frequency of. these spontaneously occurring 
almost-sinusoidal oscillations agreed very well with the value predict- 
ed by the Tollmien-'-Schlichting theory. Their amplitude increased down- 
stream until turbulent motion ensued, at first intermittently and then 
continuously. Since the occurrence of these regular oscillations was 
unexpected, numerous tests were made to rule out the possibility that 
they were effects of vibration or the result of disturbances traveling 
upstream from the turbulent part of the boundary layer or the result of 
acoustic phenomena. All these causes were ruled out. 

The reason why these phenomena had not previously been observed 
was made clear by increasing the stream turbulence tenfold to 0.2 per- 
cent. The oscillations were found, but if the turbulence had been much 
larger they would have been difficult to identify because of the 
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presence of irregular fluctuations and the near coincidence of their 
point of appearance with the transition point. In all known previous 
experiments the stream turbulence was of the order of 0.5 to 1«0 per- 
cent. Under these conditions transition was controlled by the magni- 
tude and scale of the free stream turbulence in accordance with Taylor i s 
theory. 


Systematic Experimental Study at National Bureau of Standards 

The next step was to produce waves in the laminar layer under con- 
trolled conditions and study their behavior. I will pass over the many 
schemes tried before completely satisfactory results were, obtained with 
only a mention of sound, both pure notes and random noise, from a loud 
speaker in the tunnel wall and from a email headphone behind a hole in. 
the plate. These attempts are described in the original NACA advance 
confidential report .issued in April 19^3, now declassified and available 
from the NACA as a War Time Report. Excitation by intense sound waves 
does produce transition nearer the leading edge of the plate, and, in the 
lowest turbulence conditions in the wind tunnel the fluctuations were 
shown to be mainly wind tunnel noise. In the free 'atmosphere it may 
well be that the principal initial disturbances are sound waves from the 
engine or propeller. . ' 

Schubauer and Skrams tad finally devised an ingenious method of intro- 
ducing small disturbances of known frequency as illustrated in figure 3* 
The disturbances were produced by vibrating a very thin (0.002 inch) 1 
flat ribbon whose mean position was about 0.006 inch from the: plate. 

A magnetic field was produced in the vicinity of the ribbon by electro- 
magnets on the opposite side of the plate, and an alternating current 
of the desired frequency passed through the ribbon. The amplitude of 
the disturbance could be ;cbntrolled by varying the current. The en- 
suing fluctuations of speed in the boundary layer were measured by a 
hot wire anemometer. Frequency, damping, and wave length could be 
independently measured. 

The effect of the ribbon on the mean flow was exceedingly small and 
could not be detected 2 inches downstream. The results did not depend 
on the exact dimensions or position of the ribbon. The usual procedure 
was to measure u' at various distances for a fixed wind speed and : 
frequency, then repeat at several different frequencies. From such a 
series of curves the neutral curve showing frequencies of waves which 
are neither damped or amplified could be determined and the damping 
and amplification coefficients could be determined for the. region of 
Reynolds numbers covered. The wave lengths were determined by connect- 
ing the input to the ribbon to one' pair of plates of a cathode-ray 
oscillograph and the hot-wire output to the other. The Q di stance through 
which the hot wire was moved to change the phase by 180 is one-half the 
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length 611 ^' 7116 Velocity 18 e< l ual to frequency times the wave 


Figure 4 shows the experimental values of the wave lengths of dis- 
turbances which are neither damped nor amplified corresponding to the 
Reynolds numbers indicated. Disturbances lying with.in the open loop 
are amplified; those outside it are damped. For Reynolds numbers below 
disturbances of all wave lengths are damped* 


Lin’s Revision of Theory 

the experimental work had been completed, C. C. Lin, then 
a Chinese student at the California Institute of Technology and now 
professor at Brown University, undertook a revision of the mathematical 
theory of the stability of two-dimensional parallel flows and a clarifi- 
cation of some features of the Tollmien-Schlichting theory which had 
een adversely criticized. Lin made a brief report to the National 
Academy of Sciences (reference 14) and later published a complete ac- 
count (reference 15 ). Lin’s calculation gave slightly different re- 
sults from those of Schlichting. Both curves are shown in figure 4. 

The agreement between theory and experiment is remarkable and similar 
agreement is found for the independently determined frequency locus 
J.iig. 5J. Further checks were made of the distribution of the u- 
f luctuations across the boundary layer. These showed the node and the 
phase reversal in the outer part of the layer as computed by Schlichting 
quantitative agreement. The validity and applicabilitv 
o e eories of Tollmien, Schlichting, and Lin are beyond question. 


Confirmation of Experiments at Guggenheim Aeronautical 
Laboratory, California Institute of Technology 

Early progress reports on this work were made available to H. W. 
Liepmann at the California Institute of Technology, who was engaged 
in the study of transition on the convex and concave sides of a curved 
plate. Liepmann was able to observe and study the Tollmien-Schlichting 
^waves on the convex side of the plate. Schlichting had computed the 
+ v i • conVex cur va ture on the neutral curve and found a slightly 
ev v llZinS effect ' Liepmann used the oscillating ribbon technique of 
Schubauer and Skramstad as well as methods using an artificial regular 
roughness of variable wave length, and acoustic excitation. The re- 
sults were found to agree closely with those obtained at the National 
Bureau of Standards, the influence of convex curvature being extremely 
small. The results were originally issued as NACA advance confidential 
report ACR No. 3H30 in August 1943 (now declassified and available as 
an NACA War Time Report) , 
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Relation of Tollmien-Schlichting Waves to Transition 

The mechanism of the instability of the laminar boundary layer can 
be said to be fully understood. Whatever small disturbances are ini- 
tially present are selectively amplified until large sinusoidal oscilla- 
tions are present. The small disturbances may arise from noise, free— 
stream turbulence, or surface roughness. The large oscillations do not, 
however, constitute turbulent flow. 

Some experiments were made to study the passage from waves to tur- 
bulent flow as described in the report of Schubauer and Skrarastad. The 
regular waves grow in amplitude, then become very distorted, then bursts 
of high frequency fluctuations occur. In records at some distance from 
the surface the bursts appear in the low velocity part of the cycle and 
in some records very near the surface there is evidence of intermittent 
separation. 

The general belief is that the small scale irregular eddies of the 
turbulent motion arise from some dynamic instability like the breaking 
of water waves or the rolling-up of a vortex sheet. These dynamically 
unstable vortex sheets may arise from intermittent separation during a 
part of the wave cycle, or from modification of the velocity distribu- 
tion by the shearing stress existing because of the u*v* correlation 
in the Tollmien-Schlichting waves. These stresses increase with in- 
creasing amplitude of the wave. 

It seems clear that further theoretical and experimental research 
is needed on this difficult non-linear problem. 


Effect of Free— Stream Turbulence 

The effect of free -stream turbulence at low turbulence levels has 
already been mentioned and illustrated in figure 1. Liepmann made a 
few measurements of the effect of turbulence levels from 0.0006 to 
0.003 (|j- values, v* and w’ not measured) on transition on the con- 

cave side of a curved plate. Otherwise, no work in this field has come 
to my attention. 


Effect of Curvature 

Liepmann completed an extensive study of the effect of curvature 
on transition. The work on a convex surface has already been mentioned 
as well as the result — the boundary layer on convex surfaces, at least 
up to values of momentum thickness equal to 0.001 times the radius of 
curvature, exhibits the same Tollmien-Schlichting instability as the 
flat plate and the effect of curvature iB so small as to lie within the 
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precision of measurement. Transition on the upper surface of an airfoil 
will therefore not he noticeably influenced by the curvature. 

The experimental work on concave surfaces was described by Liepmann 
in NACA advance confidential report ACR No. kJ2& (now declassified and 
available as an NACA War Time Report), When the momentum thickness 9 
is greater than 0.0005 times the radius of curvature r . the laminar flow 
is dynamically unstable due to three-dimensional disturbances as studied 
theoretically by Gortler (reference lb). In a stream of the lowest turbu- 
lence the G'6rtler parameter R 0 Jd/r, where R 0 iB the Reynolds number 
based on momentum thickness, is equal to 9*0. This value decreases to 
about 6.0 at a turbulence level of 0,003 (value of u*/u). 

If the value of 9/ r is less than 0,00005 the flow is unstable be- 
cause of Tollmien-Schlichting. waves like the flat plate and the convex 
surface and transition occurs at a constant Reynolds number whose Value 
depends on the free stream turbulence. For values of 9/ r between 
0.00005 and 0.0005 there appears to be a more or less continuous change 
from transition due to Gortler vortices to transition caused by Tollmien- 
Schlichting waves. i 


Effect of Pressure Gradient 

Schubauer and Skrametad made a few measurements of the effect of 
pressure gradient. For moderate gradients, that is, 2 or 3 percent of 
the velocity pressure per foot acting on boundary layers less than 1/2 
inch in thickness the favorable gradient had no noticeable effect on the 
neutral curve whereas the adverse gradient expanded the region of ampli- 
fication. For gradients about five times as great, comparable to those 
which might be found on airfoils, the favorable gradient gave damping up 
to the highest Reynolds number reached (about 2600 based on displacement 
thickness) and the adverse gradient gave amplification over the region 
studied, which was limited by the onset of transition. These results are 
in agreement with the effects commonly observed on transition, namely, 
that a falling pressure delays or prevents transition whereas a rising 
pressure brings about early transition. This is the phenomenon used in 
the design of laminar flow airfoils. 

Liepmann made some measurements <?f the effect of pressure gradient 
on the critical Reynolds numbers for the boundary layer on the convex 
surface of a plate with 20— foot radius of curvature. These measurements 
showed changes of the critical Reynolds number for this plate of about 
20 percent for quite moderate pressure gradients, adverse gradients (de- 
celerated flow) reducing and favorable gradients (accelerated flow) in- 
creasing the critical Reynolds number. The rate of change was most 
rapid near zero pressure gradient. For the concave surface of a plate 
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with 2.5 foot radius of curvature the effect of pressure gradient appear- 
ed to he smaller than the errors of measurement which were for this case 
rather large . 


Effect of Suction Applied to Boundary Layer 

: J. Ackeret (reference 17) and his colleagues at the Aerodynamic 
Institute of the Federal Institute of Technology at Zurich have demon- 
strated that the transition from laminar to turbulent flow in the bound- 
ary layer can be prevented by removing a part of the boundary layer air 
by suction through a small number of slots. He was in fact able to re- 
store a turbulent boundary layer to the l ami nar condition. By careful 
attention to the detailed design of the internal ducting of the suction 
system he has obtained extremely low drag coefficients. 

Mr. ,E. F. Relf (reference 18) gives a general account of work car- 
ried out in England during the war, on boundary layer suction as applied 
to the maintenance .of laminar flow and. -reduction of drag of rather thick 
airfoils. A single slot was usedon each surface of an airfoil especial- 
ly designed to suit the suction, an idea due to A. A. Griffith and fol- 
lowed up by S. Goldstein and by Lighthill on the theoretical side and 
presumably by Re If and his collaborators on the experimental side. 
Reference is also made to the use of distributed suction acting through a 
porous surface with theoretical work by J. H. Preston and a flight test 
by F. G. Miles. Wo experimental results are given and no references to 
detailed reports . 

. The information now available on the aeronautical research work 
carried out in Germany during the war (reference 19) shows that the sub- 
ject of boundary layer suction as a means of preserving laminar flow re- 
ceived attention. Various theoretical studies were made of the stabil- 
ity of 1 an y! par flow with pressure drop and with homogeneous suction. 

I have not had the opportunity of examining any of the reports in detail. 
The theoretical calculations indicate that the stability limit can be in- 
creased from its value of about 420 without suction to about 70,000 for 
the asymptotic velocity distribution finally attained at a great distance 
from the leading edge with uniform suction. 

It is hoped that details of all of this work may soon be generally 
available. 
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MECHANICS OF THE TURBULENT BOUNDARY LAYER 
Statue at the 1938 Congress 


While the 1938 Congress demonstrated a considerable advance in the 
theoretical and experimental aspects of isotropic turbulence and the ef- 
fect of free stream turbulence on transition, there was no corresponding 
advance in the theory of fully developed turbulent motion .involving shear 
flow. At the 193^ Congress Von Herman had reviewed the turbulence prob- 
lem and summarized the various physical concepts underlying the several 
theories of turbulent interchange, including Reynolds apparent stresses, 
the. Boussinesq eddy viscocity, the Prandtl momentum transfer theory, the 
Taylor vorticity transfer theory, and the Von Karman theory. Mention was 
made of Burgers* attempt to apply the principles of statistical mechanics 
and of the measurements of the turbulent fluctuations by Wattendorf and 
Reichardt in a straight channel. In the interval between 1931 and 1938 
Taylor developed his statistical theory of turbulence which was so fruit- 
ful in treating the problem of isotropic turbulence. Von Karman had 
made a generalization of this theory and attempted to treat the problem 
of shear flow. Otherwise , there were no developments which seem at this 
time to be of outstanding importance. 

At the 1938 Congress itself C. B. Millikan gave a critical discus- 
sion of various formulae proposed for the velocity distribution and 
friction in channels and circular tubes attempting to bring them into a 
consistent system. I described a method devised by H. K» Stramstad for 
direct measurement pf the turbulent shearing stress by means of a hot 
wire anemometer. Von Karman considered the problem of parallel shear 
motion from the viewpoint of ^ the theory of correlation coefficients. 

In the. course of his paper he discussed the measurements of Wattendorf 
and Reichardt previously menti one d and remarked that over the middle 
region of half of the channel u J2 was app roxim ate ly proportional to 
the shearing stress which in turn equals - u*v s . This might, he said, 
indicate the existence of a kind of statistical similarity in this region. 
Von Karman also pointed out the difficulty of reconciling Reichardt -s 
measurements of the plane stress tensor with the picture that a scalar 
mixing length determines the transfer of momentum. In the discussions 
Tollmien and Prandtl introduced the idea that the turbulent fluctuations 
might consist of two components, one derivable from a harmonic function 
and the other satisfying an equation of the heat- conduction type, that 
is, a non-diffusive and a diffusive component or a viscosity— independent 
and a viscosity dependent type. These were the only contributions bear- 
ing on the mechanics of the turbulent boundary layer. 
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Development Arising from the Methods of Buri and Gruschwitz 


Practically all of the methods available for the computation of the 
flow in turbulent boundary layers are of the "one dimensional" type in- 
volving the application of the Von Kerman integral equation for the mo- 
mentum and the use of empirical relations obtained from a few experimental 
studies of flow in convergent and divergent passages. The earliest pro- 
cedure was that suggested by Buri (reference 20). He assumed that the 
boundary layer at any section could be characterized by the momentum thick- 


ness 9 - 


n 6 ~ 

= / ^1 — an< l fcy a single parameter r fixing 


the 


shape of the velocity distribution curve. The Von Karman momentum equa- 
tion was then applied to give an ordinary differential equation connect- 


ing 9 and r • Buri first suggested a parameter 




d(T/Ty) 

d(y/0) 


where t 


is the shearing stress au distance y and T w is the shearing stress 

at the wall, that is, T a is the non-dimensional slope of the shearing 

stress curve. From the equations of motion — = e- = - pUU* at the 

dy ox 

puu f e 

wall. Hence r a = — t ♦ The shearing stress at the wall was assumed 

1 w 

to be of the form 


t w /i /2 pu 2 = f(ue/ v , r a ) 


But for T a = 0, f ~ (U0/u)~ 1/1,4 . Then Buri assumed that in the general 
case f = i (r a ) (U0/u) 1 ^ 4 , that is, that the variables are separated. 
With this empirical assumption 


-pUU»0 

r 

a i/ 2 P u 2 e(r a )(ue/v) 1/4 

It is then more convenient to use P = — l/ZTat (P a ) = 


V , 9/W\ 1/4 
[ — ) as the 

U v v ' 


form parameter. The introduction of these assumptions in the momentum 
relation gives the differential equation 


ik'r-Kc -(?*■)'} 
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In which the symbol H represents the ratio of the displacement, thick- 

6 

ness 5* = / ^1 — ^^dy to 0. On Buri*s assumptions the right-hand side 

! -'o 

is a function of T only. From his own experiments on accelerating 
flows and from Nikuradse’e experiments on both accelerating and deceler- 
ating flows, it appeared that the right-hand side was roughly a linear 
function and the equation was integrated on that assumption. Hovarth 
(reference 2l) discusses Buri f s assumptions and gives an alternative 
method of solution. Separation occurs when r is -0.07 according to 
Buri, -0.06 according to Howarth, and between -0.05 and -0.09 according 
to Prandtl, 


At about the same time that Buri proposed the method just described, 
Gruschwitz (reference 22) also suggested that the non-dimensional velocity 
distribution curves, u/U versus y/0, belong to a one-parameter family. 

He selected as the parameter q = 1 — ( — ) . Then 

v U / y=0 

„6/e - 

f = l : * s 

0 

should be a function of t). From experimental results, Gruschwitz found 
empirically for q < 0.8 

0 d 2 

: — x — (l/2pU q) = o.00894ti - 0.00461 

l/2pU 2 dx 


The momentum equation gives the relation 


T w 

pU 2 


d0 0U» 
dx U 


(2 + H) 


If some assumption is made for T w , the equations can be solved for 0 end 
q • Gruschwitz made the very rough assumption that t w is constant and 
suggested that separation takes place at or soon after q = 0.8. 

Stuper (reference 23) found "fair agreement with Gruschwitz f s equa- 
tions for observations in free flight under conditions where no separa- 
tion occurred, but Peters (reference 24) found greater and greater depar- 
tures as separation was approached. 

Von Doerihoff and N. Tetervin (reference 25) gave a comprehensive sum- 
mary and analysis of experimental data and suggested a'modified procedure. 
They selected H as an appropriate shape parameter, plotting the values 
of u/U for a series of constant values of y/0 against H as abscissa. 
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A compilation of data from many sources gave a reasonably close approxi- 
mation to a single family of curves. The value of H for separation 
vas found to lie between 1.8 and 2.6 as compared with Gruschwitz s s 
value of I.85. 


Von Doenhoff and Tetervin then proposed to use the parameter 

2fflU 9 / ^ their notation — — — — ^ as characteristic of the external 
Tv \ 1 to T o / 

forces acting on the boundary layer. This parameter is the same except 
for the sign and the factor 2 as the r a originally proposed by Buri. 

The parameter is the ratio of the nonr-dimensional pressure gradient to 
the shearing stress coefficient at the wall. It is then assumed that 
this parameter determines not H itself, that is, the particular veloc- 
ity distribution curve of the one-parameter family, but the rate of 


„dH , , .. , , . 2pUU*0 

~ is plotted against — > 


w 


1 w 


being de- 


change of H. Thus 
t6rmined from the Squire and Young formula (reference 26) 

— = |" 5.890 log lQ ( 4.075 —-)] 

The experimental data are fairly well represented by the empirical 
equation 


dH 

dx 


4 . 680 (H— 2 . 975 ) L gpUUt0 


I 


'w 


-2.035(H-1.286) 


The two equations just given taken with the momentum relation give, on 
elimination of v w twc simultaneous first order differential equations 
which can be solved by a step-to-step calculation. 

Kalikhman (reference 27) describes a different empirical approach 
based again on the momentum relation and the Squire and Young formula. 

He introduces a variable 



which reduces the momentum equation to an equation linear in z. The 
variation of the coefficients is said to be small in a turbulent boundary 
layer and the solution of the resulting linear equation with constant 
coefficients is given by the formula 


z = 
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The physical intrepretation of the assumptions made appears to he 
that the pressure gradient has a considerable influence on the resist- 
ance law and on the parameter H but only a negligible influence on U0/y. 

Kalikhman suggests that separation occurs when 

. , v-0.08 

=0.0013, 

U \ u / 

this simple criterion being an empirical approximation to a much more 
complicated expression. 

These are the principal developments along this path which have 
come to my attention. Although of some importance to the practical 
engineer, they do not satisfy the inquiring mind attempting to construct 
a rational theory. All rest on more or less arbitrary empirical assump- 
tions, and become more and more unsatisfactory as separation is ap- 
proached. 


The Fediaevsky Method 

Fediaevsky (reference 28) called attention to the fact that the pre- 
vious treatments had assumed either a constant shearing stress or a lin- 
ear variation of shearing stress through the boundary layer. These dis- 
tributions were known to be incorrect when a pressure gradient is present* 
He proposed the polynomial distribution given by 

— = 1 + A i (y/6) + A 2 (y/5) + A (y/&) + A 4 (y/&) 

Tw 


where the coefficients A are determined from the relation at the wall 
^ 3 T 

ix = F* 8114 relations at the outer boundary t = 0, Bt /B y = 0, 

B 2 T/By 3 = 0. They are found to be 

* 8 Bp . 6 Bn 

1 * t w a*' 2 = ' = " 4-3 ^ a? ’ : 

Fediaevsky proceeds by using the Prandtl assumption 


8 Bp 

T w 


= p I s ( & 

V d.v. 


A 2 


dy/ 

The mixing length l is approximated by the relation 

- = 0.14 — 0.08 (l - y/6) 2 - 0.06 (l - y/8) 
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an empirical formula developed by Prandtl on the basis of experiments 
on flow in pipes. The solution of the resulting equation in explicit 
form is a complicated equation of many terms. 


Developments Proceeding from Reynolds Theory of Turbulent Stresses 

It is not necessary to review here the fundamental equations of tur- 
bulent flow developed by Reynolds which differ from the Navier Stokes 
equation only in the additional stress components 


~ Pu |2 , - f>v e2 , - pv 1 -PU f v^ -pv^w^ and — pu f w*. 

The great success of the Taylor— von Karman statistical theory as applied 
to isotropic turbulence leads to confidence in its fundamental validity. 
Isotropic turbulence can, however, be described very simply by two param- 
eters, an intensity parameter and a scale parameter. In non— isotropic 
turbulence the description of the state of the turbulence becomes much 
more complex. Six quantities instead of one are required to specify the 
intensity and six scalar functions ^are required to specify the correla- 
tion tensor. As noted by von Karman at the last Congress the single 
available experimental measurement of a turbulent stress tensor was in- 
consistent with the idea of a simple mixing process dependent on one 
scalar quantity. 

A recent paper by Nevzgljadov (reference 29) suggests a type of 
ifheory intermediate between a complete theory utilizing a stress tensor 
and a mixing length tensor anjl the well— known simple mixing length 
theories of Prandtl and von Karman. Nevzgljadov writes down the 
Reynolds equations, the continuity equation, and the energy equation for 
the turbulent velocity fluctuations. The unknown functions are the mean 
velocity, the mean pressure, the turbulent pressure, it =1/3 p (u ,2 +v |2 4w $2 ) 

the turbulent stress tenBor, the current density of the kinetic energy 
of the turbulent fluctuations, the energy flux transferred by the pres- 
sure fluctuations, and the viscous dissipation. It is proposed to select 
the mean velocity, the mean pressure, and the turbulent pressure as in- 
dependent fundamental quantities and to assume that the other unknowns 
can be related to these fundamental quantities by "equations of Btate". 

This procedure is a generalization of the assumption made by Prandtl 
as to the relation between the . turbulent shearing stress and the mean 
motion, a relation which corresponds to what Nevzgljadov describes as an 
equation of state. The difference is that in Nevzgljadov ! s theory the 
turbulent shearing stress may depend on the turbulent pressure as well as 
on the mean motion. In fact, the question is discussed as to whether 
the shearing stress at a given point can be expressed in terms of the 
values of the fundamental quantities iiTthe immediate neighborhood of the 
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point. It might depend on their values upstream or even at all points in 
the field and thus be capable of representation only by an integral equa- 
tion. 


The particular "equation of state" connecting the turbulent shearing 
stress and the independent quantities proposed by Nevzgljadov is 

/ dui + ^k\ 

= - "( sj Si ;) 

where e, is a numerical constant. This amounts to assuming that the 
Prandtl mixing length is proportional to the turbulent velocity fluctu- 
ation. I will report at once that our experiments to be described in 
the next section do not support this assumption and our main interest in 
the paper at this time is its general philosophical background. These 
same experiments show that the turbulent shearing stress depends much 
more on the local turbulent pressure than on the local mean velocity 
gradient. 

National Bureau of Standards Experimental Program 

The experiments referred to are part of a long range 6tudy of the 
mechanics of a separating turbulent boundary layer conducted by the 
National Bureau of Standards in cooperation with, and with the financial 
support of, the National Advisory Committee for Aeronautics. Emphasis 
is placed on the measurement of the turbulent stress tensor by the use of 
the hot-wire anemometer. The problem is a difficult one because the size 
of the hot-wire configurations cannot be reduced below about 2 or 3 milli- 
meters and hence a very thick boundary layer is required. Furthermore, 
the curvature of the surface must be sufficiently small to avoid appreci- 
able pressure difference across the thick boundary layer. After several 
unsatisfactory experimental arrangements the desired conditions were 
obtained on a partition of airfoil— like section in the 10— foot wind 
tunnel of the National Bureau of Standards. This partition was 27.9 
feet long and 2 fpet thick at its maximum thickness extending in a diam- 
etral plane across the wind tunnel. The leading edge has a radius of 
1 inch . and is joined tangentially to cylindrical surfaces of 23-foot 
radius forming the nose. The trailing portion of the partition is flat 
on one side and has the form of a circular cylinder of 3^— foot radius 
on the other. This unsymmetrical shape was found necessary to secure 
separation. A blister was added at the tunnel wall to modify further the 
pressure distribution to cause separation to occur well upstream from the 
trailing edge. By various tricks two-dimensional flow was obtained over 
the central region with separation occurring uniformly at 25.7 feet from 
the nose . 
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Figure 6 shows the variation of the velocity Just outside the “bound- 
ary layer along the length of the plate. Because of the unsymmetrical 
shape of the plate at the trailing edge, the flow around the nose is very 
turbulent and produces early transition. The pressure distribution re- 
sembles somewhat that on an airfoil. These and all other measurements 
were made at a speed of about l6l feet per second at the 17— foot position, 
the Reynolds number being held as constant as possible from day to day. 

The boundary layer was about 2j inches thick at the 17-foot position and 
about 9 inches thick at separation. The thickness of the boundary layer 
at 17 feet was equal to that which would have prevailed on a flat plate 
11.3 feet long with fully turbulent boundary layer and no pressure gra- 
dient, the equivalent flat plate Reynolds number for thiB position being 
14,300,000. The turbulence in the free stream of the tunnel is about 0.5 
percent. 


A contour map of the mean speed distribution is shown in figure 7 
and the displacement thickness, momentum thickness, and their ratio, 
the parameter H, are plotted in figure 8. The velocity U in figure 7 
varies with x as shown in figure 6. The value of H at the separation 
point is approximately 2.7- Figure 9 shows the distribution of mean speed 
plotted in the manner suggested by von Doenhoff and Tetervin. To avoid 
confusion no comparison with their results is shown on the figure. The 
agreement is in general good, but there are some systematic differences 
slightly greater than the experimental dispersion, which are perhaps to 
be ascribed to differences in Reynolds number.' and pressure gradient or to 
other dif ferences in the experimental conditions. 


The National Bureau of Standards experimental program contemplates 
an exhaustive study of the turbulent fluctuations within the boundary 
layer whose general characteristics have Just been described. The first 
attention was given to the direct measurement of turbulent shearing 
stress by the method described at the last Congress. A detailed dis- 
cussion of the technique is outside the scope of this paper. In brief, 
the hot-wire head for shearing stress measurement consists of two tung- 
sten wires about 0.00031 inch in diameter and !§■ to 2 millimeters long, 
set approximately 90° to each other and each making an angle of 45° to 
the flow direction. The instantaneous voltage on one wire is proportional 
to Au + Bv; on the other to Au - Bv. Hence the root mean square volt- 
age a of the first wire is proportional to A 2 u ,a + 2AB uTv* + B 2 v’ 2 j 
for the second wire the root mean square voltage b is A 2 u |2 - 2AB uv + 
B 2 v' 2 . Hence a - b = 4AB uv and the shearing stress T is found from 


T 


= - p uv = - 


P (a - b) 
4AB 


The root mean square voltages are measured by suitably compensated 
amplifiers . 
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The measured distributions at 22^ feet and 25 feet are shown in fig- 
ures 10 and 11, in which the ordinate is the coefficient of shearing 
s tree 8 defined by 

C T = T/(l/2)(pU 2 ) 

The curves given by Fediaevsky’s polynomial approximation are also shown. 
For the 17-foot station (curves not shown), where there is no pressure 
gradient, the agreement is fairly good, and for the 25~ foot station, 
where the shearing stress at the wall is nearly zero, the agreement is 
excellent. At intermediate stations tbe agreement is poor. 

These curves give a fair idea of the scatter of the experimental data 
which are perhaps not yet so accurate as one might desire. The measure- 
ments are very tedious. It may be noted that the viscous shearing stresses 
are wholly negligible in the regions studied because of the large Reynolds 
number. The laminar sublayer is extremely thin and is never approached 
in any of the measurements . At the 17-foot position at 0.1 inch from 
the surface the turbulent shearing stress is about 190 times the viscous 
shearing stress. 

The distribution of the turbulent shearing stress has been determined 
at some 15 sections between the l4-foot station and a station just beyond 
the separation point at 25.7 feet, A contour map of C T is shown in 

figure 12. The contours show some irregularities whose causes are not 
definitely known. Beginning at the region of minimum pressure (maximum 
velocity outside the boundary layer), the maximum value of the shearing 
stress is found at greater and greater distances from the wall until at 
separation it is located near the middle of the boundary layer. The value 
of the maximum stress increases as the distance along the wall increases, 
although the actual stress does not increase anywhere near as much as the 
stress coefficient which is defined in terms of the maximum velocity just 
outside the boundary layer at the station in question. This velocity de- 
creases by a factor of about 0.7 in passing from the 17- foot to the 25- 
foot station (see fig. 6), and hence a stress coefficient based on the 
velocity at 17 feet as a constant reference velocity would be equal for 
the 25-foot station to approximately one-half the values plotted in fig- 
ure 12 for that station. The contours plotted on this basis are shown 
in figure 13. The contours are thus those of equal absolute values of 
the stress. Comparison of this diagram with figure 7 shows some sim- 
ilarity between the two. The contour for u/U =0.8 is almost identical 
with the outer contour for C Tq = 0.003; the contour for u/U = 0.6 

lies approximately along the locus of maximum shearing stress. 

The estimated value of the shearing stress coefficient at the wall 
is shown in figure l4. These values were obtained from curves like those 
shown in figures 10 and 11 by extrapolation, use being made of the fact that 
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at the wall 8 t /by = 8p/dx. The Squire— Young formula is not a good rep- 
resentation of this curve. 

The distribution of the turbulent shearing stress coefficient for 
the 17§-. 20- , and 22|-foot stations is shown in figure 15. 

The Prandtl mixing length l can he computed from the relation 


l = 

1 du 
U dy 

Typical values for the 17h> 20, and 22^-foot stations are shown in fig- 
ure l6. These curves indicate that l is a definite function of the 
distance from the wall fot y/0 up to about 1.0. 

Figure 17 shows the distribution of mean speed at these three stations 
determined by pitot tube measurements. 

Figure 18 shows the stress coefficient plotted against the mean 
velocity gradient. These curves certainly do not suggest a very defi- 
nite relation between the turbulent shearing stress and the mean velocity 
gradient . Over extended regions the velocity gradient is nearly constant; 
yet the shearing stress changes by a factor of 4 or 5. 

Surveys of u* , v* , and w* are in progress. Values for the 17i~> 

20-7, and 22^-foot stations are shown in figures 19, 20, and 21. The 
turbulence is three-dimensional, even though the mean flow is two-dimen- 
sional. The turbulence is strongly non— isotropic, v* being considerably 
less than u’, and w' being intermediate in value but closer to u*. 
Isotropy is reached in the free stream. 

In order to test Nevzgljadov’s assumption that the shearing stress 

t is proportional to (u* 2 + v* 2 + w* 2 ) — ■ , the ratio of t to u ,s 

dy 

+ v' 2 + w !2 was plotted against The results are shown in figure 22. 

j.. 

Far from being proportional to — , the ratio is almost independent of 

— except in the outer part of the layer, where T falls to zero and 
dy 

the turbulence decreases to that of the free stream. Perhaps if the free 
.'Stream were sufficiently free from turbulence, both quantities would fall 
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to zero together, retaining a constant ratio. It is evident that 
Nevzglja4ov‘B assumption is not a good approximation, hut confirmation is 
given to Von Karmen *s assumption of statistical similarity between the 
fluctuations at different points in the same cross section. The ratio 

of t to u ,2 + v ,2 + w* 2 is more nearly constant than its ratio to 
u’ 2 + v* 2 or to u* 2 or to v* 2 . 

From the data given, it is possible to compute the direction of the 
principal axis of the turbulent stress tensor by the relation 


tan 2 a 


2 u’v 


u» - v»‘ 


C T IT 


u 


t* _ v t c 


where a is the angle between the principal axis and the direction of 
flow. The results are plotted in figure 23 along with Reichardt’s re- 
sult for two-dimensional flow under pressure between two plates as given 
by Von Karman. It is seen that the results agree in general with those 
of Reichardt, The variation across the boundary layer is of the order 
of 5° to 10 °. These results emphasize the difficulty pointed out by 
Von Karman, namely, that the principal axis of dilatation is at 45° to 
the mean flow as compared with 10° to 30° for the principal axis of the 
turbulent shearing stress. Thus there are directions for which there is 
a shearing stress but no rate of shear of the mean flow and vice versa. 

I suggest that the data are to be interpreted as meaning that the turbu- 
lent shearing stress is not determined by the local mean velocity and its 
derivatives but by the local non— isotropic turbulence, and that this local 
turbulence cannot be dependent solely on the local mean velocity and its 
derivatives. The turbulence must depend on conditions upstream. 

The process of turbulent separation is seen not to depend on processes 
within the laminar sublayer, for the turbulent stress near the wall falls 
regularly and smoothly as separation is approached. At these high Reynolds 
numbers the laminar sublayer serves only as a medium for transferring the 
stress from the fluid to the wall. A further study of the conditions 
near the separation point must await the completion of the u’, v', and 
w' measurements in this region. The difficulties are great as the u' 
component becomes a largo fraction of the mean speed, and intermittent low- 
frequency disturbances are much more frequent. 


This is the point to which our present study of the data has brought 
us. The conclusions so far drawn are mostly negative in character; 
Fediaevsky 's approximation is not entirely satisfactory, and Nevzgljadov's 
assumptions are not supported by the experimental data. The Prandtl end 
Von Karman formulas relating the turbulent shearing stress to the local 
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mean motion may 1)0 fundamentally in error not only in the difficulty of 
reconciling the experimental data with the concept of a scalar mixing 
length varying from point to point but also in its fundamental philosophy. 
The stress may not be determined by the local mean motion; the local non- 
isotropic turbulence, which is directly related to the shear stress, is 
itself not determined by the local mean motion but by the generation of 
turbulence upstream, and its. convection, diffusion, and dissipation on 
the path to the point under study. 

There is, however, considerable support for Von Karman's concept of 
the statistical similarity of the turbulent fluctuations. 

The objective of the National Bureau of Standards long range experi- 
mental program is to obtain sufficient detailed information on the turbu- 
lence itself to permit the test of various theoretical assumptions and to 
suggest others which have an experimental foundation. When the u', v*, 
and w’ measurements are completed., the correlation tensor of second 
order will be studied. Unfortunately, we have not yet developed methods 
of measuring a large number of quantities at a large number of points in 
the flow field simultaneously . The measurements require months to com- 
plete. They are therefore not suitable for the analysis of differential 
changes from point to point with high accuracy. They are intended rather 
to give an over-all survey of relations between quantities directly measured 
with moderate accuracy . It is hoped that this progress report will prove 
suggestive to other workers in this important field. 


National Bureau of Standards, 

Washington, D. C., September 1946. 
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Figs. 1,3 
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Figure 1.- Reynolds number of transition as a function of free stream tur- 
bulence for a flat plate in a flow without pressure gradient. 
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Fig. 3 
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Figure 3.- Experimental arrangement of Schubauer and 

Skramstad for exciting controlled disturbances in 
the laminar boundary layer of a plate by means of a 
vibrating ribbon. 
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Figure .4.- Neutral curve (instability boundary) obtained from observed wave 

lengths for zero damping or amplification compared with the theoretical 
curves of Schlichting and Lin. Wave length » 2n/ai 6* = displacement thickness of 
boundary layer; R = Reynolds number of boundary layer based on 6*. 
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Figure 5.“ Neutral curve (instability boundary) obtained from ob- 
served frequencies for zero damping or amplification 
compared with the theoretical curve of Schlichting. 
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Figure 7,- Contour map of mean speed distribution in tbe 

turbulent boundary layer studied at NBS. u * local 
speed at point x,yj U » speed just outside the boundary layer 
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Figs. 8 % 9 




Figure 9.- Variation of u/U with H for various values of y/0. 
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Figure 8.- Displacement thickness 6*, momentum thickness 0, 
and parameter H * 6*/0 for NB8 turbulent boundary 

layer. 
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Figs. 12,13 



Figure 12.- contour map of distribution of turbulent shearing stress coefficient for 

NBS turbulent boundary layer. C T = t/1/2pU 2 ; t= turbulent shearing stress 
at point x,y; U = speed just outside boundary layer varying with x as shown in figure 6. 



Figure 13.- Contour map of distribution of turbulent shearing stress for NBS turbulent 
boundary layer. C r 0 = r/l/2pu2* u q * constant reference speed (speed just 
outside boundary layer at 17 1/2 ft station). 



















Figure 18.- Relation between shearing stress coefficient 
and the mean velocity gradient for 17 1 / 2 , 20 
and 22 l/2 ft stations. 





















